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We present a numerical calculation of the weak localization peak in the magnetoconductance 
for a stroboscopic model of a chaotic quantum dot. The magnitude of the peak is close to the 
universal prediction of random-matrix theory. The width depends on the classical dynamics, but 
this dependence can be accounted for by a single parameter: the level curvature around zero magnetic 
field of the closed system. 
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I. INTRODUCTION 

Random-matrix theory (RMT) makes system- 
independent ("universal") predictions about quan- 
tum mechanical systems with a chaotic classical 
dynamicaiSii^. The presence or absence of time-reversal 
symmetry (TRS) identifies two universality classes. 
RMT is also capable of describing the crossover between 
the universality classes, e.g. when TRS is broken by 
the application of a magnetic field B. The crossover is 
predicted to depend on a single system-specific param- 
eter, being the mean absolute curvature of the energy 
levels Ei around B = Q. More precisely, a universal 
magnetic-field dependence of spectral correlations is 
predicted when B is rescaled by the characteristic field 



dB^ 



B=0 



-1/2 



(1.1) 



with A the mean level spacing. This prediction has been 
tested in a variety of computer simulation a^'^iT . 

In open systems there exists a similar prediction of uni- 
versality for transport properties, but now the character- 
istic field depends also on the conductance g of the point 
contacts that couple the chaotic quantum dot to electron 
reservoirgSiSiifliii. A universal magnetic field dependence 
is predicted if B is rescaled by B^^, provided g is large 
compared to the conductance quantum /h. To provide 
a numerical test of this prediction is the purpose of this 
paper. 

We present a computer simulation of the open quan- 
tum kicked rotatorii^Jiiii, which is a stroboscopic model 
of a quantum dot coupled to electron reservoirs by ballis- 
tic point contacts. The ensemble averaged conductance 
increases upon breaking of TRS, as a manifestation of 
weak localization. The height, width, and lineshape of 
the weak localization peak are compared with the pre- 
dictions of RMT. 

The simulation itself is straightforward, but the for- 
mulation of the model is not. There exist several ways 
to break TRS in the closed kicked rotator^^d^dTd? and 
related modelsi^iSLSiiS^. When opening up the system 
one needs to ensure that the scattering matrix satisfies 



the reciprocity relation 

S{-B) = S'^{B) 



(1.2) 



(The superscript T indicates the transpose of the scat- 
tering matrix S.) We also require that TRS is broken 
already at the level of the classical dynamics (as it is in 
a quantum dot in a uniform magnetic field). Finally, we 
need to relate the TRS-breaking parameter in the strobo- 
scopic formulation to the flux enclosed by the quantum 
dot. All these issues are addressed in Sees. II and III be- 
fore we proceed to the actual simulation in Sec. IV. We 
conclude in Sec. V. 



II. TIME-REVERSAL-SYMMETRY BREAKING 
IN THE OPEN KICKED ROTATOR 

A. Formulation of the model 

The kicked rotator is a particle moving along a circle, 
kicked periodically at time intervals ro^i^. The strobo- 
scopic time evolution of a wave function is given by the 
Floquet operator !F. In addition to the stroboscopic time 
tq and the moment of inertia /, which we set to unity, F 
depends on the kicking strength K and the TRS-breaking 
parameter 7. We require 



(2.1) 



which guarantees the reciprocity relation (|1.2() for the 
scattering matrix when we open up the model. 

We will consider two different representations of J- ^ 
both of which can be written as an M x M unitary ma- 
trix. The classical limit corresponds to a map defined on 
a toroidal phase space. The difference between the two 
representations is whether TRS breaking persists in the 
classical limit or not. The simplest representation of J- 
has one kick per period. It breaks TRS quantum me- 
chanically, but not classically. This would correspond to 
a quantum dot that encloses a flux tube, but in which 
the magnetic field vanishes. A more realistic model has 
TRS breaking both at the quantum mechanical and at 
the classical level. We have found that we then need a 
minimum of three kicks per period. 
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B. Three-kick representation 



is then given by 



We will mainly consider the three-kick model, so we 
describe it first. In this model TRS is broken both quan- 
tum mechanically and classically. Stroboscopic mod- 
els with multiple kicks per period of different magni- 
tude were studied previously in the context of quantum 
rachets^ii. Inspired by that work, we study the time- 
dependent Hamiltonian 



H{t) = P- + ^V{e)J2Mt^n + e) + 5,{t-n-e)] 

11 

+ 7Cos(6i)^(5(t-?i + l/3) 

n 

-7Cos(6i)^(5(t-?i- 1/3), (2.2) 

n 

with e an infinitesimal. The angular momentum operator 
p = —ificsde is canonically conjugate to the angle E 
[0, 27r). The effective Planck constant is fieff = ^tq/I. 
The potentialiLiSiSiSi 

V{0) = K cos(7rg/2) cos(6l) -I- sin(7rg/2) sin(26l) (2.3) 

with q 7^ breaks the parity symmetry of the model. The 
form of the potential is such that in the large X-limit the 
diffusion constant does not depend on q. For 7 = there 
are two kicks per period in Eq. (|2.2I) . but since they are 
displaced by an infinitesimal amount we still call this a 
"single-kick" model. For 7^0 two more kicks appear 
with opposite sign at finite displacement. We will see 
that this choice guarantees the reciprocity criterion (|2.1f) 
for the Floquet operator. 

The reduction of the Floquet operator 



T = T exp 



^off Jo 



H{t)dt 



(2.4) 



(with T the time ordering operator) to a discrete, finite 
form is obtained only for special values of ?ieff, known 
as resonances^^. We have to reconsider the usual condi- 
tion for resonances in the presence of additional, TRS- 
breaking kicks. Here our analysis departs from the quan- 
tum rachet analogySS. 

The initial wave function ?/'(^) evolves in one period to 
a final wave function given by 

^{B) = exp(-zl/(0)/2?ieff)exp(zfieff5|/6) 

X exp(— z7C0s(6')/?ieff) exp(z?ioffC^e/6) 

X exp(i7Cos(6')/?icff) exp(i?icff9|/6) 

X exp(-i\/(6l)/2neff)V'(6'). (2.5) 

One recognizes three factors describing free propagation 
for 1/3 of a period, each followed by a kick. The res- 
onance condition for free propagation is ?ieff = 27rr/M, 
with T an odd integer and M an even integer—. The free 
propagation 



27r 
3m' 
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E 

m,n'— 



27rr r 
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3M^ V V 3Af 



(2.7) 



Resonance means that the initial and final wave functions 
can be treated as discrete vectors on a 3Af-point lattice, 
labeled by the indices n, n' . The angle 9 is an arbitrary 
offset parameter. Different values of are not coupled by 
the free propagation. Putting together three iterations of 
Eq. H2.7|l we get three independent components of + 
27rn/3Af) for n = 0,1,2 (mod 3), each on an A/-point 
lattice. 

We find that the resonance property is preserved in 
the presence of intervening TRS-breaking kicks, provided 
that r = 3 and M even, but not a multiple of 3. The 
free propagation ()2.7|l then is conveniently expressed in 
matrix notation. The matrix acts on an M-component 
vector Vm = i'iO + 2Trm/3M), m = 0, . . . , M - 1. We 
choose the arbitrary phase ^? = 0, so that 



mm' Yn 

rn'— 

The matrices are defined by 



A —irrm jM 

j^,j~l/2^-27Timm' /M 



(2.8) 



(2.9) 
(2.10) 



The matrix product J7^E[/ can be evaluated in closed 
form, with the result 

= M-i/2e'"/4exp[i(7r/M)(m' - m)% 

(2.11) 

Collecting results, we find that for hcs — 6n/M the 
Floquet operator (|2.5I) is represented by an M x M uni- 
tary matrix, of the form 



Y 

rr, 



{xnY*iiYnx) 



„i(M7/67r) cos(27rm/M) 
-j(M/127r)y(27rm/M) 



(2.12a) 
(2.12b) 
(2.12c) 



V'i(0) = exp(i?ieff5gVg)^(^) 



(2.6) 



One readily verifies the reciprocity relation (|2.1|) . 

The classical map corresponding to this quantum me- 
chanical model is derived in App. A. We show there that 
TRS-breaking of the classical map is broken for 7 7^ in 
the three-kick model. 



C. One-kick representation 

TRS breaking in the one-kick model is constructed 
as a formal analogy to the magnetic vector potential. 
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by adding an offset 5 to the momentum of the kicked 

rotatorl5-16.17.18^23,24^^ 

To obey reciprocity 

T{-S)=T^{S) (2.13) 

for odd M it is enough to symmetrize the expression of 
Ref. 15 by infinitesimaUy sphtting the kick (as it was 
done in Ref. jL4 for S — 0). For even M, which is more 
convenient for appHcation of the fast Fourier transform, 
one also needs to redefine the lattice points in order to 
preserve reciprocity^Si. 

The model takes the form 



Urr. 



= M 



1/2 -27rj(m~ 



I'/M 



-i(MK/A-iT) cos(27rm/M 
M 



T(m— - 



-1 ,M, 



/M 



(2.14a) 
(2.14b) 
(2.14c) 

(2.14d) 



In addition to the TRS-breaking phase 5 there is a phase 
<f> to break the parity symmetry. The reciprocity property 
(|2.13ll can easily be checked. 

The classical map corresponding to this model is also 
discussed in App. A. It does not break TRS. 



D. Scattering matrix 

To model a pair of iV-mode ballistic point contacts that 
couple the quantum dot to electron reservoirs, we im- 
pose open boundary conditions in a subspace of Hilbert 
space represented by the indices ml"''. The subscript 
n = 1,2,...N labels the modes and the superscript 
a — 1,2 labels the leads. A 2A^ x M projection ma- 
trix P describes the coupling to the ballistic leads. Its 
elements are 



P„. 



1 if rn = n e {rrir"^}, 
otherwise. 



(2.15) 



The mean dwell time is tjj = M/2N (in units of Tq). 

The matrices P and together determine the scatter- 
ing matrij 



.12.13.14 



S{e) = P[e-'' - Til - P^P)]-^TP^, 



(2.16) 



where e is the quasi-energy. The reciprocity condition 
(|2.1|) of implies that also S satisfies the reciprocity 
condition (|1.2|l . 

By grouping together the N indices belonging to the 
same point contact, the 2A^ x 27V matrix S can be de- 
composed into 4 sub-blocks containing the N x N trans- 
mission and reflection matrices, 



r t 
t' r' 



(2.17) 



The conductance G (in units of /h, disregarding spin 
degeneracy) follows from the Landauer formula 

G = Trtt^ (2.18) 



III. RELATION WITH RANDOM-MATRIX 
THEORY 

In RMT time-reversal symmetry is broken by means 
of the Pandey-Mehta Hamiltonian^i 



H = Ho + iaA, 



(3.1) 



which consists of the sum of a real symmetric matrix 
Hq and a real antisymmetric matrix A with imaginary 
weight ia. We denote by Mh the dimensionality of the 
Hamiltonian matrix. The two matrices Hq and A are 
independently distributed with the same Gaussian dis- 
tribution. The variance i/^ = ((-ffo)?,) — {^fj) (* 7^ j) 
determines the mean level spacing A = -nv j ^JMu at the 
center of the spectrum for Mh 3> 1 and a ^ 1. 

To lowest order in perturbation theory the energy lev- 
els Ei{a) depend on the TRS-breaking parameter a ac- 
cording to 



SE, 



42 



E, 



(3.2) 



with SEi = E,{a) - Ei{0) and E, = Ei{0). The char- 
acteristic value Qfc is determined by the mean absolute 
curvature, 



-1/2 



(3.3) 



From Eq. I3.2|l we deduce that ac — A/z/ ~ 1/\/Mh, 
up to a numerical coefficient of order unity. A numerical 
calculation gives 



'Mh = krmt = 1-27. 



(3.4) 



A real magnetic field B is related to the parameter a 
of RMT by 



(3.5) 



where Be is determined by the level curvature according 
to Eq. (|1.1() . For a ballistic two-dimensional billiard (area 
A, Fermi velocity vp) with a chaotic classical dynamics, 
one has^i^ 



B, = c—iAVA/hvFy/^ 
eA 



(3.6) 



with c a numerical coefficient that depends only on the 
shape of the billiard. The field Be corresponds to a flux 
through the quantum dot of order [h/ e)^TcrgA/Ti ^ 
h/e, with the ergodic time being the time it takes 
an electron to explore the available phase space in the 
quantum dot. 

The analogue of Eqs. (|1.1|) and (|3.5|l for the quantum 
kicked rotator considered here is 




^72 



-1/2 



(3.7) 
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FIG. 1: The critical value 7c of the TRS-breaking parameter 
in the closed three-kick model is presented for different system 
sizes at fixed K. The parity-breaking parameter is g = 0.2. 
The solid line shows the large- A" limit 13.811 . The dashed lines 
are averages over M of the numerical data. 
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FIG. 2: Same as Fig. Q but now for the closed one-kick 
model. The parity-breaking parameter is = 0.2 tt. The 
solid line shows the large- A' limit ^A.'s)\ . 



Here 7 is the TRS-breaking parameter in the three-kick 
model. The same relation applies to the one-kick model, 
with 7, 7c replaced by 5, 5c- 

To complete the correspondence between the kicked 
rotator, RMT, and the real quantum dot, we need to 
determine the two characteristic values 7c and 5c- In 
App. B we present an analytical calculation deep in the 
chaotic regime [K — s- 00), according to which 

lim 7c = 127rM-3/2KRMT = 47.9 Af "^/^ (3.8) 
lim 5c = 4V3Af"3/^KRMT = 8.80 AT^/'^ - (3.9) 

In Figs. n and 121 we show a numerical calculation for finite 
K , which confirms these analytical large- -fC limits. 

In the open system the characteristic field scale for 
TRS-breaking is increased by a factor with g the 
conductance of the point contacts. We consider ballistic 
iV-mode point contacts, so that g — N, measured in rmits 
of e^/h. The conductance G{B) of the quantum dot is 
also measured in units of e^/h. According to RMT, the 



GiB) 



N 
1 



1 r 



1 + (2KRMTiV-l/2B/Bc)2 



.9.11 



(3.10) 



For the quantum kicked rotator we would therefore ex- 
pect a weak localization peak in the conductance given 
by 



G(7) = G„ 



1 + (2KRMTiV"'/'7/7c)' 



(3.11) 



in the three-kick model. We define the weak localization 
correction 5G{^) — G{j) — Goo, with Goo the conduc- 
tance at fully broken TRS. The expression in the one-kick 
model is similar, with 7/7C replaced by d/5c- 

In the large- limit we can use the analytical expres- 
sions (|3.8|) and H3.9|) for 7c and 5c to obtain 



lim 5G{j) = -i[l + (M3/2iV-i/2^/67r)2]-i, (3.12) 
+ (Af3/2Ar-i/2^/2V3)2]-i.(3.13) 



K 

lim dG{5) 

K^oo 



In App. C we show how these two results are consistent 
with a semiclassical calculation. 



IV. NUMERICAL RESULTS 

The numerical technique we use to calculate the con- 
ductance was described in Refs. 14 and |2^ The calcu- 
lation of the scattering matrix (|2.16(l is performed ef- 
ficiently by use of an iterative procedure and the fast- 
Fourier-transform algorithm. We need to average over 
many system realizations (varying lead positions and 
quasi-energies) to suppress statistical fiuctuations. In ad- 
dition, we need several points to plot the 7-dependence. 
This makes the calculation for large M more time con- 
suming than earlier studies of universal conductance fiuc- 
tuations in the same model at zero magnetic ficld^^'^^. 

First we present in Figs. 3 and 4 results for the weak 
localization correction 5G in the three-kick model as a 
function of the TRS-breaking parameter 7. The data 
arc obtained by averaging over 40 lead positions and 80 
quasi-energies. The parameter 7c was calculated for the 
closed model using Eq. H3.7|) , and the resulting RMT pre- 
diction H3.11|) is also shown (dotted curve). 

To compare the simulation with RMT in more detail 
we have fitted a Lorentzian 



5G 



-i[l + (A//7/7*)'i 



(4.1) 



to each data set. This is the RMT result (ICT) if 7* 
7^MT = 7c^^/^/(2^/2TB'KRMT). The large if-hmit is 



lim 7rmt = 6vr/ 



(4.2) 



In Fig. [SI we plot the fitted crossover parameter 7* as a 
function of M for fixed dwell time. The plot confirms 



the scaling with 



1/2 



cx g- 



and also shows good 
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FIG. 3: Dependence of the average conductance on the TRS- 
breaking parameter 7. The three-kick model is characterized 
hy K ^ 7.5, q = 0.2, and r p = M /2N = 25. The dotted line 
shows the RMT prediction 13.1111 . with 7c calculated from the 
mean level curvatures (Fig.0. 



FIG. 5: Dependence of the crossover parameter 7* on the 
system size. The data are obtained by fitting the Lorentzian 
14.111 to the numerical data of Figs. |3| and 2] The solid line 
shows the large A'- limit 13.81 . The dotted lines are the RMT 
prediction for K = 7.5 and K = 41, using 7c found from the 
level curvatures in the closed model (Fig. . 
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FIG. 4: Same as Fig. El but for a: = 41. 



agreement with the values of 7rmt calculated from the 
mean level curvature (dotted lines). 

We also performed numerical calculations for the one- 
kick model. The crossover scale S* extracted from a 
Lorentzian fit to the weak-localization peak was com- 
pared with the value (5J^mt = ^c-M^^^/(2\/2tdkrmt) pre- 
dicted by the mean level curvature. The large iC-limit of 
this value is 



lim S 

K- 



RMT 



= \/6/ 



(4.3) 



We show in Fig. O the ratio ^*/<5^mx for the one-kick 
model, as well as the ratio 7*/7rmt f*-"^ ^^'^ three-kick 
model. The ratio is close to unity for both models if the 
dwell time is sufficiently large. At the smallest to there 
is some deviation from unity in the one-kick model. 

The magnitude of the weak localization peak in Figs. 
13 and 01 shows a small (about 10%) discrepancy with the 
RMT prediction. We attribute this to non-ergodic, short- 
time trajectories. We show in Fig. jSl the dependence of 
the magnitude of the weak localization peak (5G(0) on the 
dwell time. The results suggest that 5G{Q) + \ 1/td: 
a deviation from RMT to be expected from the Thouless 
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FIG. 6: Dependence of the ratio 7*/7rmt for the three-kick 
model and the ratio (5*/^r,mt for the one-kick model on the 
dwell time td ■ Data points for a given dwell time are obtained 
by averaging over system sizes in the range from 200 to 1000. 



energy scale (which is oc 1/td). The deviation from miity 
in Fig. has presumably the same origin. 

We could determine the M-dependence of 7* and 
(5G(0) up to M = 10" (for K = 7.5 and td = 5). The 
motivation for extending the calculation to large system 
sizes is to search for effects of the Ehrenfest time22*2i. 
Alt houg h the Ehrenfest time te ~ 3.8 (estimating as in 
R.ef. UM was comparable to re = 5, we did not find any 
systematic M-dependence in 7* or (5G(0), cf. Figs.[31and 



IZI 



V. CONCLUSIONS 

In conclusion, we have studied time-reversal symmetry 
breaking in quantum chaos through its effect on weak lo- 
calization. We have found an overall good agreement 
between the universal predictions of random-matrix the- 
ory and the results for a specific quantum mechanical 
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FIG. 7: Dependence of 5G{0) on the system size Ai for several 
dwell times. Dashed lines show averages over system size. 
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FIG. 8: Dependence of the amplitude of the weak localization 
peak SG{0) (averaged over several system sizes) on the dwell 
time td- Dashed lines show a linear dependence on 1/td, 
extrapolated to the RMT value 1<5G(0)| = 1/4. 



model of a chaotic quantum dot. In particular, the scal- 
ing oc g^^^^ of the crossover magnetic field with the point 
contact conductance g is confirmed over a broad param- 
eter range. 

Deviations from RMT that we have observed scale in- 
versely proportional with the mean dwell time td , consis- 
tent with an explanation in terms of non-crgodic short- 
time trajectories. These deviations therefore have a clas- 
sical origin. 

More interesting deviations of a quantum mechanical 
origin have been predicted^Si^ in relation with the finite 
Ehrenfest time te- This is the time scale on which a wave 
packet of minimal initial dimension spreads to cover the 
entire available phase space. The theoretical prediction is 
that the weak localization peak 5G{0) cx e""^^/"^" should 
decay exponentially once te exceeds T]j. Our simula- 
tion extends up to te — to, but shows no sign of this 
predicted decay. This is consistent with the explanation 
advanced by Jacquod and Sukhorukov^^ for the insen- 
sitivity of universal conductance fluctuations to a finite 
Ehrenfest time. As pointed out in Ref. the same ex- 
planation also implies that weak localization should not 
depend on the relative magnitude of te and to- 



Because our simulation could not be extended to the 
regime te > to, this final conclusion remains tentative. 
It might be that the exponential suppression of 5G{0) 
does exist, but that our system was simply too small to 
see it. 
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APPENDIX A: CLASSICAL MAP 

Here we derive the classical map that is associated 
with the quantum mechanical Floquet operator of the 
kicked rotator with broken TRS. We consider separately 
the three-kick and one-kick representation. 

1. Three-kick representation 

We seek the classical limit of the Floquet operator 
(|2.12|) . We consider the classical motion from 6*0 at i = 
to Ot a.t t = T (in units of tq). Intermediate values of the 
coordinate are denoted hy 9t, t = 0,1, ...,T. The classical 
action S is the sum 



T-l 



(Al) 



s = J2s{0t+i,et). 

t=Q 

Following the general method of Ref. ? we derive 

s{e', e) = Sc{0', 02) + Sb{02,Oi) + Sa{0i,0), (A2) 

Sa{0i,e) = 1(01 - + 2TTag, f - 67rap,0i 

+ 7cos(0i)- (A3) 
Sb{02,0i) = 1(02 - 01 + ^TTae. f ^ 6TTap,02, (A4) 
Sci0',02) = 1(0' - 02 + 27:ag,)^ -67:ap,0' 

--fcos{02)- ^V{0). (A5) 

The integers ag , are the winding numbers of a classical 
trajectory on a torus with e [0, 27r) andp G [0, Gtt). The 
map equations are derived from 

PI = ^Sa{0Y,B), p^- — Sa{0l,0), (A6) 



P2 
P' 



001 
d 

802 

d_ 

80 



Sb{02,01 
Sc{0',02), P2 



00 

Pi = -^^Sb{02,0i), (A7) 



d0' 



■Sc{0',02). (A8) 
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Eqs. ljA6IA8p are equivalent to the following set of 6 equa- 
tions that map initial coordinates {0,p) onto final coor- 
dinates {9',p') after one period: 



6*1 =e+p/3-V'{e)/6^2Trag,, 
Pi = p - 7 sin 01 - V'{9)/2 - Girap, , 

6*2 = 01 +pi/3- 2TTag^, 

6*' = 02 + P2/3 4- 7 sin 6*2/3 - 2TTa0, , 
p' ^P2+l sin 6*2 - V'{0')/2 - QnGp. . 



(A9) 
(AlO) 
(All) 



We denote V — dV/dO. Winding numbers of a trajectory 
on the torus in phase space {0,p) are denoted by ug, <Tp. 
These integers are determined by the requirement that 
9,61,62,6' e [0,27r) and p,pi,p2,p' e [0,67r). TRS for 
a classical map means that the point {6',^p') maps to 
{6,—p). This property is satisfied for 7 = 0, but not 
for 7 7^ 0. TRS is broken at the classical level in the 
three-kick model. 



2. One-kick representation 

We now seek the classical limit of the Floquet operator 
(|2.14|l . The classical action S after one kick is 

S{6',6) = ^{9' ^ 6 + 2Traef - 27rap6' 
+ 5{6' -6 + 2nag) 

- lK[cos{6 + (/))+ cos(6i' -I- (/))]. (A12) 
The map equations are derived from 



P 



d_ 

06 



-S{9',9), p. 



d_ 



s{e',9). 



(A13) 



The mapping of initial coordinates {9,p) onto final ones 
[9' ,p') after one kick is then 



9' = 9+p+^Ksm{ 



5 — 2Trag, 



p' = p + iA'[sin(6' + (/>)+ sm{9' + </))] - 2Trap 



The canonical transformation p — S ^ p, 9 
brings the map to an equivalent form 

9' = 6 +p+\Ks\u6 -2-11(78, 
p' = p+\K[s\n6 + sin9')~2TTap. 



(A14) 



(A15) 



This form is manifestly invariant under the transforma- 
tion that maps {6' , —p') onto {6, —p) for any value of 
and 5. Hence TRS is not broken at the classical level in 
the one-kick model. 



APPENDIX B: DERIVATION OF EQS. (IsTsll AND 

In the large- -fC limit the level curvature in the kicked ro- 
tator can be related to the level curvature in the Pandey- 
Mehta Hamiltonian. This leads to the relations H3.8|) and 



(|3.9|l between the TRS breaking parameters 7 (three-kick 
model) and 5 (one-kick model), on the one hand, and the 
Pandey-Mehta parameter a, on the other hand. 

Perturbation theory for eigenphases (j)i{5^) of a unitary 
matrix T{5^) gives the series expansion 

^^{h) - 0» + Wu5l + \ E,^. iW^yf ((57)'cotan^^ 
^\y^^{hf■ (Bl) 



Here 
C/diag(e'' 



denotes an eigenphase of .?^(0) 



The Hermitian matrices 



W and V are defined by = U {-iT'^ d~^T\^=n)V\ 
V — djW\j^o- Due to reciprocity of we find Wu = 0. 
For the three-kick model 12.1211 the operators W, V are 



OTT 



UC)YUXU\ (B2) 



Af 

V = i( — )2[/xtHtyt(cntCH - HtCHC)rHAC/t, 

(B3) 

where Cmm' = Smm' cos{2TTm/M). We assume that for 
strongly chaotic systems {K ^1) the matrix elements 
Wij and Vii are random Gaussian numbers independent 
of the eigenphases. Average diagonal elements calcu- 
lated in the three kick model at 7 = are {Vu) = 
TtV/M = and (Wu) = TrW/M = 0. The variance of 
the off-diagonal elements is = TrWWyAP = 

M/(67r)2. 

For the one-kick model H2.14|l the operators W, V are 



W = UX^U^DUXU\ V 



— M, (B4) 

ZTT 



with Dmm' = 5mm' {m +1/2 — M/2 — 5M /2tt). Average 
diagonal elements at (5 = are (V^i) = TvV/M = -M/2'k 
and {Wu) = TrW/M = 0. The variance of the off- 
diagonal elements is P) = TrWWyM'^ = M/12. 

For K ::§> 1 the eigenphases 0i are distributed ran- 
domly in the circular ensemble, which is locally equiv- 
alent to the Gaussian ensembleV We expand Eq. ljBl|) 
for small eigenphases difference, compare with Eq. 1)3.2(1 
and substitute the variances of matrix elements calcu- 
lated above. For the one-kick model we drop terms with 
Vii as they are of order 1/M. We finally arrive at Eqs. 
and (pTi^l . 

The explicit formula for the Pandey-Mehta parameter 
a describing the kicked rotator at large K is 



3/2 



127r 



(B5) 



for the three-kick model. The corresponding formula for 
the one-kick model is 



(B6) 
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APPENDIX C: SEMICLASSICAL DERIVATION 
OF THE WEAK LOCALIZATION PEAK 

We present a semiclassical derivation of the weak lo- 
calization peak, adopting the method of Ref. to the 
case of the kicked rotator. The method can not be used 
to determine the amplitude SG{0), but we use it for the 
crossover scale. This serves as an independent check for 
the scaling predicted by RMT. 

The action difference in the three-kick model for a pair 
of trajectories related by TRS is calculated as follows. 
The action Sq for a trajectory with initial coordinate Oq 
and final coordinate 0t at 7 = is compared with the 
action S for a trajectory with the same initial and final 
coordinates, but at small 7. The result of linear expan- 
sion in 7 is 

AS' = 5-5o =7^[cos6'i(t)-cos6i2(i)], (CI) 
t 

where periods are numbered by t — 0,1, ...,T — 1 and 
0i{t),d2{t) denote the coordinate of the particle when 
TRS-breaking kicks are applied. 
The weak localization correction is 

6G cx (exp(2iA5'/;ieff)) , (C2) 

where the average is taken with respect to all trajectories 



connecting initial to final coordinates. Approximating 
the distribution of the phase difference A^* for a single 
step by a Gaussian, and taking the continuum limit of 
exponential dwell-time probability P{t) oc e"*/"^", we de- 
rive 

5G«[l + (M7/7T]-\ h*r^2nlJ{rn^), (C3) 

with v being the variance of AS/j for a single step. The 
result = 1 for large K (and large td) is obtained by 
averaging over random initial points in the whole phase 
space. We thus find Eq. (|4.2|) . the same result as the one 
obtained in RMT. 

The action difference for a pair of symmetry related 
trajectories in the one-kick model is 

AS = S -Sa = Sj2Wit)-e{t) + 2TTae{t)], (C4) 
t 

to linear order in 5. This leads to 

<5G cx [1 + iMS/S*)^]-\ (<5*)2 = 2hlg/{TDiy). (C5) 

By averaging over random initial points in the whole 
phase space for large K and td we find ly = 47r^/3. Hence 
we obtain Eq. the resuh of RMT. 
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